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ABSTRACT 

Current observational evidence does not yet exclude the possibility that dark energy 
could be in the form of phantom energy. A universe consisting of a phantom con- 
stituent will be driven toward a drastic end known as the 'Big Rip' singularity where 
all the matter in the universe will be destroyed. Motivated by this possibility, other 
evolutionary scenarios have been explored by Barrow, including the phenomena which 
he called Sudden Future Singularities (SFS). In such a model it is possible to have 
a blow up of the pressure occurring at sometime in the future evolution of the uni- 
verse while the energy density would remain unaffected. The particular evolution of 
the scale factor of the universe in this model that results in a singular behaviour of 
the pressure also admits acceleration in the current era. In this paper we will present 
the results of our confrontation of one example class of SFS models with the available 
cosmological data from high redshift supernovae, baryon acoustic oscillations (BAO) 
and the cosmic microwave background (CMB). We then discuss the viability of the 
model in question as an alternative to dark energy. 
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1 INTRODUCTION 



The d is covery of the acceler ation of the universe l|Riess et al.l 
1998), (|Spergel et al.l [20031 ) resulted in the breaking of the 
link between geometry and destiny of the universe. It was 
now the nature of the substance driving this acceleration 
that would determine the ultimate fate of the universe. Cur- 
rent observations are consistent with a cosmological con- 
stant as the origin of this acceleration, but they are not yet 
able to rule out other possibilities such as phantom energy 
which will drive the universe towards a 'Big Rip' singular- 
ity. This situation has encouraged the study of other ex- 
otic singularities. One such example is the Sudden Future 
Si ngularity (SFS) (Type II acco r ding to the classification 
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the Big Separation (Type IV (Dabrowski & Dcnkicwic: 
20101) ) and the ui-singularity (jDabrowski fc Denkiewic: 
2009). These singularities may serve as alternatives to dark 
energy (|Dabrowski fc Denkiewicz! |2010| ). Furthermore, the 
Big Brake is a special kind of a singular ity of the SFS type, 
which can arise in tachyonic models (|Kamenshchik et al.l 
In this paper we present the results of our investi- 
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Noiiri. Odintsov fc Tsuiikawal (120051 )). first proposed by 
arrow ] ( 20041 ). Other exotic singularities discovered to date 



include the Finite Scale Factor (FSF) singularity (Type III 
l|Noiiri et al.l I200E )), the Generalised Sudden Future Sin - 
gularity (GSFS) {Dabrowski. Denkiewicz fc Hendrvl [20071 ) . 
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gation of one example class of models which accommodate 
an SFS, also proposed by Barrow, by confronting it with the 
current cosmological observations. 

The present work is not the first time that the SFS 
models have bee n confronted with cosm ological observa- 
tions. Previously, iDabrowski et all (|2007h performed such 
a test using the luminosity distance redshift relation ap- 
plied to supernovae data. Their paper showed that these 
data were consistent with the SFS model over the redshift 
range probed by the supernovae, and moreover could per- 
mit an SFS that would occur in a suprisingly short time: 
less than 10 million years from now. In this paper we extend 
the analysis to confront this class of SFS models with some 
other available cosmological observations and constraints, 
from the cosmic microwave background radiation, baryon 
acoustic oscillations and the age of the universe. 

The structure of this paper is as follows. In Section 2 we 
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introduce some relevant theory underlying Sudden Future 
Singularity models. Section 3 gives an account of the cos- 
mological probes we have used to constrain our SFS model, 
followed by an explanation of the methods we employed to 
perform these tests. In Section 4 we present the results of our 
investigations. Finally in Section 5 we give our conclusions. 



2 SUDDEN FUTURE SINGULARITY MODELS 



In 2004 iBarrowl 1)20041 ) first published the results of his dis- 
covery of the existence of Sudden Future Singularities, which 
arise in the expanding phase of a standard Friedmann uni- 
verse and violate the dominant energy condition only. This 
was an intriguing discovery since until then the theoretical 
search for expanding universes with possible violent ends 
had identified only Bi g Rip singular ities which violate all 
the energy conditions ||B arrow! |2004| ) . Consider the Fried- 
mann equations: 
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where p is the energy density, p the pressure, a the scale 
factor, k the curvature index, c the speed of light, G the 
gravitational constant and an overdot represents derivative 
with respect to time. Barrow found that by keeping the scale 
factor and the Hubble expansion rate finite one will also nec- 
essarily maintain a finite density, but that the pressure can 
still diverge. The divergence of pressure is accompanied by 
the blow up of the accelera tion, leading to a scala r polyno- 
mial curvature singul arity ([Ha wking & Elli slll973lV Shortly 
after their discovery, iNoiiri fc Odintsovl (|2004al lbl) showed 
that these singularities may be avoided (or moderated) when 
quantum effects are taken into account. 

The finiteness of the energy density and the fact that 
p — y oo at an SFS means that the null (pc 2 + p ^ 0), weak 
(pc 2 ^ and pc 2 +p ^ ), strong (pc 2 +p and pc 2 +3p > 
0) energy conditions are satisfied but the dominant energy 
condition (pc 2 ^ 0, ~pc 2 ^ p ^ pc 2 ) is violated. 

We can see that the divergent behaviour of the pressure 
cannot be linked to the finite energy density or we will not 
have an SFS. Hence we need to release the assumption of 
an equation of state which imposes a link between these two 
quantities. It should be noted that an SF S can occur also i n 
inhomogeneou s and anisotropic models <|Dabrowskill2005h , 
iBarrowfooolh . 

Similar to sudden singularities are the so called 'qui- 
escent' singularities w hich occur in braneworld models 
l|Shtanov fc Sahml 120021 '). These are a type of sudden sin- 
gularities whereby the pressure remains finite alongside the 
density and the Hubble p arameter while higher derivatives 
of the scale factor diverge. lAlam fe Sahnil l|2006t ) confronted 
braneworld models with quiescent singularities with SNe la 
and BAO data and concluded that they would not fit obser- 
vations. Furthermore the exact same sudden future singular- 
ity with the div ergent pr e ssure occurs in nonlocal cosmol- 
ogy as found bv lKoivistol (|2008T ). The accelerating universe 
in this model will lead it towards an SFS rather than a de 



Sitter type epoch. Koivisto however shows that these sin- 
gularities may be avoided by a slight modification of the 
nonlocal model. 

Physically sudden future singularities manifest them- 
selves as momentarily infinite peaks of tidal forces, but 
geodesic completeness is satisfied in our SFS model. Hence 
Sudden Future Singularities are regarded a s weak singu- 
larities (|Fernandez-Jambrina fc Lazkoj [20041 ) which means 
in- falling ob servers or detectors would not b e destroyed by 
tidal forces (|Dabrowski fc Denkiewiczl 1201(1 ). The universe 
will continue its expanding evolution beyond such weak sin- 
gularities until for example the occurrence of a more serious 
singularity that is geodesica lly incomplete like the Big Rip 
which can end the universe (|Dabrowski et al]|2007l ). 

Barrow constructs an example SFS model where the 
scale factor takes the form: 



a(t) = A + Bt m + C(t s - t) r ' 



(3) 



where A > 0, B > 0, m > 0, C and n > are free constants 
to be determined. By fixing the zero of time, i.e. setting a(0) 
= 0, and the time of the singularity, a(t s ) = a s , the scale 
factor may be written in the equivalent form: 



a(t) = A + (a s - A)(^) r 
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iDabrowski et al] (|2007h changed the original parametrisa- 
tion for the scale factor by using A = Sa s , to obtain: 



i(t) = a 3 [5+(l-5)y m -S(l-y) n ] 



V = T , (5) 



where a s , n, m, S and t s are constants to be determined. This 
way they created a non-standardicity parameter, S, which, 
as it tends to zero, recovers the standard Friedmann limit 
(i.e. a model without an SFS). 

For a pressure derivative singularity to occur we need: 



1 < n < r 



r — integer. 



(6) 



We note that, for r ^ 3, all energy conditions are ful- 
filled. These singularities a re called Generalised S udden Fu- 
ture Singularities (GSFS) (|Dabrowski et al.|[2007l ) and may 
occur in theories with higher-order curvature corrections 
l|Noiiri et alj|200j ). For an SFS we need r = 2 which means 
that 1 < n < 2. 

An important requirement is that the asymptotic be- 
haviour of the scale factor close to the Big Bang singularity 
follows a simple power law clbb = y m which will simulate 
the behaviour of flat k = barotropic fluid models with 
m — 2/3(w + 1). This will ensure that all the standard ob- 
served characteristics of the early universe such as the CMB, 
density perturbations and Big Bang nucleosynthesis are pre- 
served. All the energy conditions are satisfied if we require 
m to lie in the range 

< m < 1. (7) 

Furthermore, in accordance with (7), if for a standard dust- 
dominated universe we take m = 2/3 our SFS model will re- 
duce to the Einstein-de-Sitter (EdS) universe at early times. 

The other parameters of the model that we should con- 
sider are: a s , t s and 8. The parameter a s , which sets the 
physical size of the universe at the time of the SFS, will 
cancel out in the equations of the standard cosmological 
probes we have used to test the model. Thus we do not need 
to consider this parameter further. 
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For constraining the time, t S) when an SFS might occur, 
we can introduce the dimensionless parameter yo = to/t s , 
where to is the current age of the universe. Since the sin- 
gularity is assumed to be in the future, it follows that 
< y < 1. 

The current acceleration obtained in the SFS model, 
as a result of the particular form adopted for the scale fac- 
tor, leads to a divergence of the pressure at sometime in 
the fu ture evolution of the model. iDabrowski fc Denkiewic3 
refer to the cause of this late time acceleration in the 
SFS model as a 'pressure-driven dark energy'. 

With the parameters n, m and yo having definite val- 
ues or ranges, it remained to identify a suitable range 
of investigation f or th e non-standardicity parameter 8. In 
IDabrowski et all l|2007l ) negative values of S were associated 
with acceleration, but we revisited this question in order 
to check rigorously the range of values of 8 which should 
be considered, taking into account all relevant physical con- 
straints. In so doing we imposed the established observa- 
tional facts that both the first and second derivatives of 
the scale factor are currently positive - i.e. we rejected any 
combination of SFS model parameters for which do ^ or 
'do ^ 0. In addition we imposed the physical condition that 
positive and negative redshifts should correspond to past 
and future events respectively. Finally we checked the sign 
of a throughout the evolution of the SFS model and rejected 
parameter combinations that would predict contraction (i.e. 
d < 0) at any point in the interval ^ z ^ 1, on the 
(conservative) grounds that the expansion of the universe is 
securely observed from e.g. the Hubble diagram of type la 
supernovae over this range of redshifts. We did not, however, 
make any further assumptions about the expansion history 
of the universe outside of this range; in particular parameter 
combinations that would predict e.g. future contraction of 
the universe were not excluded from our analysis. 

Thus by fixing m = 2/3, as previously discussed, and 
varying n and yo over their permitted ranges, we sought to 
identify those values of 8 which were consistent with the 
above physical constraints. It quickly became apparent that 
this task was not possible analytically. Therefore we carried 
out a numerical exploration of the 3-dimensional parame- 
ter space (n, yo,8). This then told us that 8 should not be 
positive. 



3 COSMOLOGICAL CONSTRAINTS AND 
ANALYSIS METHODS 

We carried out a comparison between our SFS model and 
current cosmological observations within a Bayesian frame- 
work, computing posterior distributions for the SFS model 
parameters inferred from each of the cosmological probes 
considered. From Bayes' theorem we may write 



p(G|data) oc p(data|G)p(G), 



(8) 



where G denotes the parameter(s) of the SFS model and 
'data' denotes generically the observed data for one of the 
cosmological probes under consideration. The term on the 
left hand side of eq. (8) is the posterior distribution for G, 
i.e. it describes our inference about the model parameters 
in the light of the cosmological data with which we are con- 
fronting the SFS model. The first term on the right hand 



side is the likelihood function which gives the probability of 
obtaining the data that we actually observed, given a partic- 
ular set of model parameters. The second term on the right 
hand side is the prior probability distribution for the model 
parameters, based on e.g. theoretical considerations and/or 
previous observations. For each of our cosmological probes 
we took the likelihood function to be Gaussian in form, i.e. 



p(data|0) oc exp(- 



1 ^ 



(9) 



where the constant of proportionality could be determined 
from the normalisation conditions that the Likelihood func- 
tion and posterior probability distribution should integrate 
to unity. In eq. (9) the quantity \ 2 takes the general form: 



X ^ (^obs ^prcd) 



(10) 



where a; b s is the observed value of the relevant cosmological 
probe, Xprcd is the value predicted by our SFS model given 
the parameters Q, a is the uncertainty (due to measurement 
error and/or intrinsic variations) in the observed value and 
n is the number of data points associated with it. We now 
briefly describe each cosmological probe in more detail. 



3.1 SNe la redshift-magnitude relation 

We calculated the redshift-magnitude relation for our model 
to compare it with the observed relation as probed by SNe 
la as standard candles. Although as mentioned in the pre- 
vious section an SFS could occur regardless of the curva- 
ture of the universe, to make our calculations simpler, and 
also to be consistent with the present observational data, 
we shall consider only a flat universe which bears an SFS. 
We used 557 SNe from the Union2 dataset as compiled by 
lAmanullah et al.l |2010) which is the largest published SNe 
la sample to date. In our analysis we treated the Hubble con- 
stant as a nuisance parameter, marginalising over a range 
of values using as a prior distribution the most recent re- 
sults from the Hu bble Space Telescope (HST) Key Project 
l|Riess et al.l 120091 ') . which assume no underlying cosmology 
in measuring this parameter. We describe this marginalisa- 
tion process in more detail below. For our SNe la likelihood 
function we adopt the \ 2 quantity 
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where [i denotes distance modulus, o" p hot is the mea- 
surement uncertainty tabulated for each SN and <7j n t is 
the scatter intrinsic to every supernova because they are 
not perfect s tandard candles . We take o\ n % to be 0.15 as 
advocated bv lKowalski et all l|2008h . 

Fig. [1] shows the distance modulus as a function of 
redshift for the SFS model and the concordance model 
for their respective best-fitting parameters and how they 
compare with the best SNe dataset available. One can see 
that the SFS model fits the data very well and that it is 
almost distinguishable from the concordance model over 
the redshift range probed by the Union2 sample. 
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R{zcm-r) = (1 + ^cmb)-Da(-2:cmb), 



(13) 
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Figure 1. The predicted distance modulus plotted against red- 
shift for the SFS model and Concordance Cosmology standard 
model (CC). The SFS model prediction is shown as the solid cyan 
line and is for parameters n = 1.995, S = —0.5, yo = 0.805. (These 
are the best-fitting SFS parameters for the Union2 SNe dataset). 
The concordance model prediction is shown as the dashed pink 
line and is for p arameters Q m = 0.27 25, Qa = 0.7275 (from the 
WMAP7 results llKomatsu et alj|201lh '). Both sets of calculations 
assu me Hp = 74.2 km s^ 1 Mpc - 1 from the HST Key Project re- 
sults (Ricss ct al. 2009). Also shown are the observed data points 
(with quoted l-cr errors) for the Union2 SNe dataset. One sees 
that the SFS model fits the data very well and that the fit is al- 
most indistinguishable from that for the concordance model. (A 
colour version of this figure is available in the online journal.) 



3.2 CMB distance priors: R and l a 

In order to test a model against the cosmic microwave 
background (CMB) data one would ideally calculate the full 
angular power spectrum of the temperature anisotropies 
for the model, and compare it with the observed angular 
power spectrum. However, since this process is rather 
computationally intensive and complex, a simpler approach 
is instead to calculate the distance scales to which the 
power spectrum is very sensitive. The positions of the peaks 
and troughs of the CMB power spectrum, which can be 
measured precisely, provide a measure of the distance to 
the decoupling epoch. The distance ratios measured by the 
CMB are: 

(1) the angular diameter distance to the last scatter- 
ing surface (at redshift zcmb, which we take to be 1089), 
Da(zcmb), divided by the sound horizon size at the decou- 
pling epoch, r s (zcMB) which is quantified by the 'acoustic 
scale' and is defined by: 



la = (1 + ZCMB 



tvDa(zcmb) 
r s (zcMB) 



(12) 



where the factor (1 + zcmb) accounts for the fact that 
Da(zcmb) is the proper ang ular diameter distance and we 
calculate j",,(zcmb) following I Wang fc Mukheried (|2007h . 

(2) the angular diameter distance to the last scatter- 
ing surface divided by the Hubble horizon size at the 
decoupling epo ch, which is cal l ed th e 1 shift parameter' , R 
and is given bv lKomatsu et al l (|2009T i: 



where f2 m is the ma t ter de nsity parameter. It was shown by 
I Wang fc Mukheried (|2007T l that using these two parameters 
together is necessary in order to place tight CMB constraints 
on the parameters of the model of interest. They found that 
models with the same parameter R but different values for 
l a , and vice versa, in general gave different CMB angular 
power spectra. 

For the CMB constraints we therefore computed likeli- 
hoods with the following % 2 expressions for the shift param- 
eter and the acoustic scale respectively: 
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The observed values for these quantities: R = 1.725 ± 

0.018 and l a = 302.09 ± 0.76, are obtained from the 
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WMAP7 data and give n by [ Komatsu et al.l 11201 if). 
ever, as pointed out by lElgarov fc Multamakil |2007j) the 
shift parameter is not a directly measurable quantity and it 
is in fact derived from the CMB data assuming a specific 
cosmological model. Care therefore needs to be taken when 
using this quantity as a cosmological constraint. 

The value for t he shift parameter quoted above, as 
iKomatsu et al l (|2009h explains, has been derived assuming a 
standard FLRW universe with matter, radiation, curvature 
and dark energy components. The SFS model we are con- 
sidering is also a standard Friedmann model but it assumes 
no explicit dark energy component; instead cosmic acceler- 
ation is driven by the divergence of pressure resulting from 
the particular form of scale factor adopted in the model. 
Matter is also permitted in the SFS and in fact, since we re- 
quire our model to reduce to the EdS case at early times, we 
adopt the same matter content as that in the concordance 
model. C oncerning the radiatio n and curvature components, 
we follow IKomatsu et all {2009) and ignore these in the shift 
parameter calculation. 

Turning to the dark energy componen t , here we fol- 
lowed the approach of lElgarov fc Multamakil J2005h and ex- 
pressed our SFS model in a form equivalent to an evolving 
dark energy model by computing its effective equation of 
state, u> e ff (z). (See the Appendix for a short derivation of 
the expression for w e ff). 

We then looked at the evolution of this effective w in- 
dex to see how it compared with the observed behaviour. 
Fig. [2] shows the evolution of the effective equation of state 
over the redshift range < z < 20, for two representative 
sets of SFS model parameters. In both cases we see the the 
same general features: w e a ~ — 1 as z — > and w e a — > 
for large z. These limiting behaviours are in good agree- 
ment with current observations. Note, however, that in plot 
(b) the particular SFS model parameters result in the di- 
vergence of u> e ff at certain redshifts, which is caused by the 
denominator in the expression on the right hand side of eq. 
A. 5 tendi ng to zero at these red shifts. This behaviour is dis- 
cussed in IShafieloo et"aH l|2006h where they make the case 
that in models where dark energy is not treated as an explicit 
fluid or a field, the equation of state cannot be used as a fun- 
damental quantity and indeed an effective equation of state 
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model in question firstly the evolution of matter density per- 
turbations during the matter dominated era must be sim- 
ilar to the Concordance Cosmology case, and secondly the 
comoving distance to the horizon at the epoch of matter- 
radiation equality should scale inversely with £l m Ho- While 
these conditions are not met for the Carneiro et al. model 
of a time varying cosmological constant, they are met in our 
case. 



Figure 2. The evolution of w a g as a function of redshift in the 
SFS model, for two representative sets of model parameters. Both 
plots show generic behaviour as z — > and z — > oo which is in 
good agreement with current observations, although plot (b) also 
shows 'fake' singularities due to the paramctrisation of w e ff. 



may display unusual properties like singul arities. Very sim- 
ilar di vergent behaviour was found e.g. bv lSola fc Stefancid 
(2005) in computing an effective equation of state for their 
evolving A model; indeed those authors refer explicitly to 
each divergent feature as a 'fake singularity, which is noth- 
ing but an artefact of the EOS parametrisation'. In our case 
too the divergence of w e a is not seen as an indication of a 
fundamental physical problem with the model. Nevertheless 
the general similarity of the limiting behaviour of w e fi to 
tha t of the concordance model gives us confidence, follow- 
ing |Mgaro^^^Multamaki| (|2007f ). that it remains appropri- 
ate to use the 'observed' value of the shift parameter when 
investigating our SFS model. 



3.3 BAO distance parameter, A 

Baryon acoustic oscillations which originate from the exci- 
tation of sound waves in the early universe photon-baryon 
plasma through cosmological perturbations are useful dis- 
tance indicators at the current epoch. We can use this to 
constrain very well the following quantity, termed as the 
'distance parameter' , using the current data: 
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(16) 



, where H(z) is 



The function E(z) is defined as: E(z) 
the Hubble expansion rate. Zbao is the effective redshift of 
the galaxy sample used to measure the distance parameter 
and it is 0.35 for the SDSS galaxy sample. For the BAO 
distance parameter constraint we adopt: 
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(17) 



Here the observed value for the distance parameter has 
been taken from the late st SDSS results which is: A — 
0.469(n/0.98)-°' 35 ±0.017 l|Eisenstein et al.ll2005h . We used 
the val ue of n = 0.963 for th e spectral index from WMAP7 
results |Komatsu et al.ll201lh . 

We find that the BAO distance parameter is not im- 
mune to the model dependency issue in it s derivation ei- 
ther. T his issue is con sidered in detail in e.g. Carneiro et al.l 
i|2008l ) who, following lEisenstein fe Hul <|l99ct) . identify two 
(implicitly assumed) conditions which should be valid in or- 
der that the BAO distance parameter is applicable. For the 



3.4 Age of the universe, to 

Using the standard Friedmann equation for calculating the 
age of the universe we have: 



t = A- 



dz 



(l + z)E(z) 



(18) 



where Dh is the Hubble distance and E(z) is defined as 
before. In our likelihood function for this constraint we used 
the x 2 expression: 

prcd \ 2 

(19) 
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where tg re is the age of the universe predicted in our SFS 
model, as computed from eq. (18 ). For to hB , the observed age 
of the universe, we followed e.g. iBalbi. Bruni fc Quercelliru] 
l|2007h and used the best current estimate derived from var- 
ious astrophysical probes, including globular cluster ages, 
that have been determined without assuming a particu- 
lar cosmological model. We adopt to = I2.6I2I2 Gyr with 
an asymmetric error distribution corresponding to 95 per 
cent upper and lower co nfidence limits, as reported in 
iKrauss fc Chaboverl (|2003T ). 



3.5 Hubble Constant, Ho 

The current value of the Hubble constant is the final con- 
straint which we employed in our analysis. Although the 
Hubble constant already features indirectly in the SNe la 
and age constraints, since it enters into our calculations of 
the predicted values for those two probes, we considered it 
useful also to compare directly the observed value of Ho with 
its predicted value calculated for our SFS model. The latter 
is simply given by: 



rrpr 
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+ n5(l - yo) 
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We thus computed a likelihood function for this constraint 
using the \ 2 quantity 



2 _ (Hq 
Xh - 



H pred-)2 



(21) 



where Hg bB — 74.2 ±3.6 km s 1 Mpc 1 , i.e. the (cosmology- 
independent) HST Key project value discussed earlier. 



4 RESULTS 

In this section we present posterior distributions for our SFS 
model parameters inferred from each of the cosmological 
probes described in Section 3. For ease of presentation we 
have computed a series of 'slices' through the 2-dimensional 
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conditional distribution of the parameters n and yo at given 
values of the non-standardicity parameter 8. We adopted 
uniform priors for n and yo over the intervals 1 < n < 2 
and < yo < 1 respectively, on the theoretical grounds 
discussed in Section 2. Thus, for these uniform priors, the 
mode of the posterior distribution function, for each 'slice' 
in 8 was coincident with the maximum of the conditional 
likelihood function. 

For the case of the SNe la data, as noted in Section 3 we 
computed the posterior distribution after first marginalising 
over the Hubble constant. It follows straightforwardly from 
Bayes' theorem that: 

p(n, yo \S) = J p(n,yo\S)p(Ho)dH , (22) 

where p(Ho) represents our prior information on the Hubble 
constant. Thus for the SNe la case we computed p(n,yo\5) 
following eq. (22), adopting a Gaussian prior for Ho with 
a mean value of 74.2 and a standard deviation of 3.6, in 
accordance with the results of the HST Key Project. 

We computed posterior distributions adopting both a 
regular grid of (n, yo) values over the 2-dimensional pa- 
rameter space and a simple MCMC-based approach using 
the Metropolis algorithm. Both methods gave very similar 
results, although the latter approach is more practical 
when exploring more general parameter spaces with larger 
numbers of parameters. We will consider such cases further 
in a subsequent paper. 

Figs. [3]- [5] present contour plots showing the Bayesian 
credible regions (evaluated at the 68 per cent, 95 per cent 
and 99 per cent level) for n and yo, computed for each 
of our cosmological probes at a series of fixed values of 8. 
In each figure we show separately contours for each of the 
six cosmological probes we have considered, labelled (a) - 
(f) respectively. These probes are: (a) the SNe la redshift- 
magnitude relation; (b) the CMB shift parameter, R; (c) 
the CMB acoustic scale, l a ; (d) the BAO distance parame- 
ter, A; (e) the present-day value of the Hubble constant, Ho; 
(f) the present age of the universe, to- Fig. [3] shows results 
for 8 = —0.7. For larger (i.e. less negative) values of 8 the 
posterior distributions for most of the probes were qualita- 
tively similar; significantly, however, the predicted value of 
the shift parameter was such that no credible region for the 
SFS model parameters was found at the 99 per cent level. 
We return to this point below. 

One can clearly see from Fig. [3] that there is no part 
of the (n,yo) plane where the credible regions for the six 
cosmological probes overlap, for this value of 8. There is 
substantial overlap between the SNe la and BAO contours; 
this is not surprising given that both probes are sensitive to 
the E(z) function over a similar range of redshifts. 

The age of the universe constraint (f) is rather weak, 
reflecting the large uncertainty on the observed value, and 
only a small region of the (n, yo) plane, on the upper left 
of Fig. HHf), is excluded at the 95 per cent level; these pa- 
rameter values are already strongly excluded by all other 
probes. 

In Fig. [3][b) we have magnified the shift parameter con- 
tours for greater clarity. Note, however, that the minimum 
X 2 value for the shift parameter in this plot is already 21.35 



(corresponding to a predicted value of R = 1.641) which is 
unacceptably large - i.e. the likelihood function (and hence 
the posterior probability) for the SFS model parameters is 
everywhere vanishingly small for the shift parameter at this 
value of 8. Similar behaviour was observed for larger values 
of 8; indeed the minimum \ 2 value becomes progressively 
higher as 8 increases. Hence we do not show contour plots 
for 8 > -0.7. 

It is interesting to note from panels (b) and (c) in Fig. 
that the two CMB constraints R and l a produce credible 
regions which do not in fact show any overlap at the 99 per 
cent level. This illustrates the importance of the point made 
by Wang and Mukherjee, as noted earlier, that one should 
use both CMB probes in order to better constrain (or indeed, 
as is the case here, to reject) a given model, since these 
probes are sensitive to the model parameters in different 
ways. Moreover we also see that the credible regions for R 
and l a do not share a common overlap with any of the other 
four probes either. Hence a fit to the SFS model parameters 
is strongly excluded for this value of 8. 

Figs. 2] and [S] present credible regions for two further 
values of 8, in order to illustrate the pattern of behaviour 
exhibited by the various probes as 8 decreases. Firstly Fig. [4] 
shows results for 8 — — 1. We can see that, as before, there 
appears to be overlap between the credible regions for SNe, 
BAO distance parameter, Hubble constant and age of the 
universe; indeed the contours for these probes have shifted 
only slightly from their position in Fig. [3] However all four 
show no overlap at the 99 per cent level with either of the 
CMB probes, which in turn continue to show no overlap 
with each other. 

Finally, Fig. [5]shows the credible regions for 8 = —1.5. 
Here we can see that the contours follow more or less the 
same shapes as in the previous figures. There is still no 
overlap between R and l a , and neither CMB probe over- 
laps with any of the other probes, again at the 99 per cent 
level. Note, moreover, that while the general appearance of 
the credible regions is similar to Fig. [4] the contours for the 
acoustic scale, l a , in Fig. [SJc) have shifted slightly down- 
wards and appear to be approaching the lower limit of yo- 
As 8 is decreased further this trend continues and hence no 
fit is obtained. To emphasize that indeed the contours dis- 
play no overlap with one another, we show in Fig. [6] the 
most important contours of SNe la, CMB shift parameter, 
R, CMB acoustic scale, l a and BAO distance parameter su- 
perimposed for values of 8 previously considered in Figs. [3]- 
[5] One can see clearly in these figures that while the SNe la 
data are consistent with an SFS occurring in the very near 
future, as shown by iDabrowski et all (|2007f l. the same SFS 
parameters would not give a fit to the CMB data. 

Note that for the purpose of straightforward illustration 
the contour plots presented in Figs. [3]- [5] have been calcu- 
lated first without the imposition of the physical constraints 
discussed in Section 2. A typical example of the application 
of these conditions, and their impact on the credible regions, 
is shown in Fig. [7] for the acoustic scale, l a . One can see in 
Fig. [TJb) how the contours are abruptly cut off, in this case 
by the inclusion of the condition that the universe currently 
be accelerating. In other words, in the region enclosed by the 
contours in the right hand part of Fig.[7Ja), e.g. for yo < 0.6 
and n > 1.5, our SFS model predicts a value of l a which is in 
satisfactory agreement with observations, but for a universe 
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Figure 3. The contours show 68 per cent, 95 per cent and 99 per cent credible regions in the (n,yo) parameter space for a fixed value 
of 5 = —0.7. The plots labelled (a), (b), (c), (d), (e) and (f) correspond to the contours calculated using the probes, SNe la redshift- 
magnitude relation, CMB shift parameter, R, CMB acoustic scale, l a , BAO distance parameter, A, the Hubble constant and the age of 
the universe respectively. 




Figure 4. The contours show 68 per cent, 95 per cent and 99 per cent credible regions in the (n,yo) parameter space for a fixed value 
of 8 = — 1. The labels are as described in Fig. [3] 
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Figure 6. The panels show a superposition of the credible regions in the (n, yo) parameter space derived from the SNe la, CMB shift 
parameter, R, CMB acoustic scale, l a and BAO distance parameter, A data, for the values of 8 = —0.7 (a), 8 = —1 (b) and <5 = —1.5 
(c), as previously illustrated in Figs. [3] -\E\ respectively. Note that in each panel there is no part of the parameter space where all credible 
regions overlap. (A colour version of this figure is available in the online journal). 



which is not currently accelerating. The other cosmological 
probes such as the age of the universe and the BAO dis- 
tance parameter are similarly affected, although much less 
severely, for certain values of 8. However, since it is already 
the case that we find no significant overlap between the cred- 
ible regions over the entire parameter space without applying 
these additional physical constraints, we will not present any 
further plots that do include them. 



5 SUMMARY AND CONCLUSIONS 

In this paper we have investigated one class of Sudden Fu- 
ture Singularity models proposed by Barrow, by confronting 




1.2 1.4 1.6 1.8 1.2 1.4 1.6 1.8 



Figure 7. The impact of the imposition of the physical conditions 
discussed in Section 2 is shown here for the l a contours. The plot 
labelled (a) shows these contours not considering the physical 
conditions while in plot (b) these conditions arc included. 
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it with the currently available observational data. After in- 
troducing the theory behind the model and explaining its 
characteristics in Section 1, we reported on the cosmological 
probes we used in testing our model in Section 2, taking care 
throughout to consider thoroughly the applicability of each 
probe to the SFS model under study. We then presented 
the results of our investigations in Section 4. Specifically 
we found that the SFS model was not compatible with all 
of the current cosmological observations considered: while 
good agreement could be found with e.g. the Union2 SNe 
data and BAO distance parameter at low redshift, the same 
model parameters predicted an acoustic scale and a shift 
parameter for the CMB which were in strong disagreement 
with WMAP7 observations. 

We have presented our results as a series of conditional 
'slices' through the (n,yo) parameter space at fixed values 
of 6. We note that the rejection of the SFS model does not 
depend on our choice of 2-dimensional parameter space on 
which to conduct our analysis. The choice of the (n, yo) space 
was convenient because of the ranges of these two parame- 
ters, but the absence of any region which gave a fit to the 
model is a robust result that can be extended to the full 
3-dimensional space of (n, yo, S). In other words, by showing 
that there is no fit to the data in conditional 2-dimensional 
spaces we have therefore shown that there is no fit to the 
data in the full 3-dimensional parameter space. 

In conclusion we note that our results are for a fixed 
value of m = 2/3, which ensures that our SFS model reduces 
to the EdS universe at early times. In a forthcoming paper 
(Dabrowski, Denkiewicz, Ghodsi & Hendry in preparation) 
we will relax this assumption and consider SFS models in 
which the parameter m is allowed to vary, by confronting 
these models with the same cosmological observations which 
we have considered here. 
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APPENDIX A: DERIVATION OF WW 

In this section we present a short derivation of the expression 
for w c ff discussed in Section 3.2. A dark energy component 
with equation of state parameter w(a), will correspond to a 
density that varies with the scale factor according to: 



p(a) = p exp -3 / — [l + w(a)] 



f a da 
Ji a 



(Al) 



The standard Friedmann equation for a flat model will 
therefore take the form: 

^) 2 = ^ + (l-On)cxp (-3^ a ^[l + W (a)]).(A2) 

Re- writing the above using a — (1 + z)^ 1 we have: 

w(z)d(z) 



x =fi m (l+z) 3 +(l-fi m )exp 



H 

which means: 



+ z 



,(A3) 



l + w(z)d(z) _ 1 / E(z) 2 - n m (i + 



1 + 2 3 ln V I - <>„ 



(A4) 



where E(z) = H(z)/Ho- Now, differentiating both sides 
gives, where ' = ^: 
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